Accepted for publication in ApJ Letters 

Preprint typeset using LAT^X style cmulatcapj v. 11/10/09 



O 

(N 

> 
O 

o 



VIB RATIONALLY EXCITED HCN IN THE LUMINOUS INFRARED GALAXY NGC 4418 
Kazushi Sakamoto 1 , Susanne Aalto 2 , Aaron S. Evans 3,4 , Martina C. Wiedner 5 , and David J. Wilner 6 

Accepted for publication in ApJ Letters 

ABSTRACT 

Infrared pumping and its effect on the excitation of HCN molecules can be important when using 
rotational lines of HCN to probe dense molecular gas in galaxy nuclei. We report the first extragalactic 
detection of (sub)millimeter rotational lines of vibrationally excited HCN, in the dust-enshrouded 
nucleus of the luminous infrared galaxy NGC 4418. We estimate the excitation temperature of T v ib ~ 
230 K between the vibrational ground and excited (i>2=l) states. This excitation is most likely due to 
infrared radiation. At this high vibrational temperature the path through the v 2 —l state must have 
a strong impact on the rotational excitation in the vibrational ground level, although it may not be 
dominant for all rotational levels. Our observations also revealed nearly confusion limited lines of CO, 
HCN, HCO+, H 13 CN, HC 15 N, CS, N 2 H+, and HC 3 N at A ~ 1 mm. Their relative intensities may 
also be affected by the infrared pumping. 

Subject headings: galaxies: ISM — galaxies: active — galaxies: individual (NGC 4418) 
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1. INTRODUCTION 

Hydrogen cyanide (HCN) is a major tracer of dense 
molecular gas in space. It decays fast in its rotational 
energy ladder, 300 times faster than carbon monoxidcQ, 
because of its large dipole moment. Rotationally ex- 
cited HCN therefore suggests an excitation mechanism 
fast enough to counter the decay, such as frequent colli- 
sions with interstellar H 2 in dense molecular gas. HCN or 
any molecule with a large dipole moment, however, will 
not trace dense gas if there is another excitation mech- 
anism that is faster than the H 2 collisions and indepen- 
dent of gas density. One such excitation path is through 
a vibrationally excited state, to which molecules can 
be pu mped by infrared radiation (|Carroll fc Gold smith 
Il981h . This excitation is possible because a round trip 
of radiative transitions between the vibrational ground 
and excited states can change the rotational level J by 
2. The first vibrationally excited state of HCN is its 
bending state (u 2 =l) 1024 K above the ground (Fig. 
[1]). Radiative excitation of HCN through this state has 
been observ ed around dusty evolve d and young high- 
mass stars (jZiurvs &: Turner] [T986). Meanwhile, de- 
spite the increasing number of extragalactic HCN ob- 
servations, it has been controversial whether such ra- 
diative excitation is significant in extragalactic obser- 
vations where the volume probed is much larger and 
having enough IR flux there more difficu lt (Aalto et al 
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found significant in a 10 14 L quasar (jWeifi et alJl2007t 
iRiechers et aHl2010| ) but the quasar may be exceptional. 

In this Letter, we report our detection of rota- 
tional transitions from vibrationally excited HCN in 
the luminous infrared galaxy NGC 4418. We ana- 
lyze HCN excitation by measuring vibrational and ro- 
tational temperatures from the data. Our detection 
follows the detection of vibrationally exc it ed H C3N in 
the same galaxy by iCostagliola fc Aaltol (|2010j ). Re- 
lated HCN lines are the ?; 2 =l-<— a bsorption at 14 /xm 
in NGC 4418 and other galaxies (jLahuis et al.l 120071) 
and the direct Z-type tra nsitions of v i bratio nally ex- 
cited HCN in Arp 220 (jSalter et al.l l2008h . NGC 
4418 has a high luminosity and a compact and heav- 
ily dust s hrouded nucleus (Lib. = lO* 1 " 1 .^, size < 100 
pc = 0"6 ; iDudlev fc Wvnn-Williamsl H99W Evans et al.l 
2003). Both an active galactic nucleus and a very young 
starburst have been proposed for the luminosity source 
(|Spoon et al.ll200ll iRoussel et al.ll200l . 

2. SMA OBSERVATIONS 

We observed the nucleus of NGC 4418 (R.A. = 
12 h 26 m 54.62 s , Dec.=-00°52'39.4"; J2000) using the 
Submillimeter Array (SMA0 in 2010 March and May 
in the subcompact and compact configurations, respec- 
tively. Receivers were tuned to around 350 GHz in March 
and 270 GHz in May for J=4-3 and 3-2 transitions of 
HCN and HCO+ as well as CO(3-2). Signals were sep- 
arately recorded from upper and lower sidebands (USB 
and LSB), whose central frequencies are 12 GHz apart 
from each other. Each sideband is 4 GHz wide; we 
filled its central 30 MHz (^30 km s -1 ) gap by interpola- 
tion. Spectral resolution was 1.625 and 3.25 MHz for our 
270 and 350 GHz runs, respectively. The system gain, 
passband, and flux scale were calibrated using a nearby 
quasar 3C273 (which dimmed in May), quasars and plan- 
ets, and Mars (T b 270G Hz=190 K and T b 3 50GHz=192 K) 

8 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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and Titan (Tb 270GHz=77 K for continuum), respectively. 

We reduced the data using MIR and MIRIAD packages 
and made images with the natural weighting of visibil- 
ities at 30 km s _1 velocity resolution. The full width 
at half maximum of the synthesized beam was 5'.' 2 for 
the 350 GHz data, which had projected baselines of 6- 
25 m, and 2"6 for the 270 GHz data (15-68 m). We 
did not subtract continuum in order to later determine 
the continuum and lines in the spectra. There are two 
calibration notes. First, we took the spectral index of 
3C273 into account in our gain calibration in order to 
accurately scale the flux of our two sidebands. We used 
the spectral index a = —1.0 (for S v oc v a \ the mean 
of our March and May measurements. Second, we veri- 
fied our passband calibration by setting aside a fraction 
of 3C273 data when determining passband shapes and 
by applying the same calibration to the galaxy and the 
3C273 data. The flatness of the resulting test 3C273 
spectra confirmed our calibration. 

3. RESULTS 

Strong continuum and line emission were detected 
at the nucleus of NGC 4418. No extended emission 
was detected off the nucleus except CO (3-2) that was 
marginally resolved. Any missing flux in our data is es- 
timated to be less than our calibration error on the basis 
of the rat io of our flux to single-dish flux, 0.82 ± 0.19 for 
CO(3-2) (lYao et al.H20"ol and 1.21 ±0.19 for HCO+(4- 
3) (Aalto et al. in prep.) 

Our spectra and emission parameters in the central kpc 
are presented in Fig. [2] and Table [1] The strongest line 
is CO(3-2) followed by HCN and HCO+ in both (4-3) 
and (3-2), CS(7-6), and N 2 H+(3-2). We also identi- 
fied H 13 CN and vibrationally excited HCN next to the 
CO and HCO + lines, respectively. Although both are 
at about +400 km s _1 from the systemic velocity of the 
adjacent brighter lines, they cannot be a high-velocity 
gas because such a feature is absent in CO (2-1) (Fig. 
[3]). This HCN is excited in its bending mode («2=1) 
and the line is one of the /-doublet (Z=lf). The other 
line (^=le) is only —38 km s _1 from the ground state 
HCN line and indistinguishable. At even lower levels, 
our data show line features that we attribute to HC 15 N 
and HC 3 N. The data plausibly contain more lines, e.g., 
on both sides of HCN(3-2). Line density is high, and we 
are approaching the confusion limit. 

The continuum, which we determined from 'line free' 
channels (details are in Table [l| , has a spectral index of 
a = 2.7 ±0.5 in our data excluding the 350 GHz USB. It 
is consist ent with preyious submillimeter measure ments 
(a ~ 2.9: IDunne et^l2000l : lYang fc Phill ips 2007). The 
350 USB continuum is an outlier with a = 6.25 ± 0.84 
between the 350 GHz sidebands (cf. a = 2.40 ± 1.04 
between our 270 GHz sidebands). The 350 USB 'line 
free' channels may thus contain blended lines. The con- 
tamination would be 23 mJy (i.e., 12%) if we scale the 
LSB continuum with a = 3. This possible error is not 
included in Table [1] and it would affect weak lines more. 

4. HCN EXCITATION 

We obtained the following HCN excitation tempera- 
tures among (w 2 ',J) = (l 1 ^,4), (l 1 -^), and (0,4) states 
under the assumptions described below: T cx (l 1 ^,4; 
1 1/ ,3)=31±13 K, T ex (l 1/ ,4; 0,4)=227±18 K, and 



1^(1^,3; 0,4)=254±24 K. The first one is a rotational 
excitation temperature T rot and the second is a vibra- 
tional excitation temperature T v jb. We needed opac- 
ity correction to the HCN(4-3) flux for these excitation 
temperatures because the line is almost certainly sat- 
urated considering its low flux ratio, 3, to H 13 CN(4- 
3). We adopted 3000±1000 Jy km s" 1 as the opacity- 
corrected HCN(4-3) flux on the basis of H 13 CN(4- 
3) and HC 15 N(4-3) fluxes and the abundance ratios 
[ 12 C/ 13 C]«50 and [ 14 N/ 15 N] «300 (|Wilsonl Il999l Milky 
Way inner disk values) . We also assumed that the vibra- 
tionally excited lines are optically thin, because if they 
were thick they would be as bright as the i>2=0 lines. An- 
other assumption here and in the T ex calculation above 
is that all HCN lines arise from the same volume. 

The vibrational excitation is most likely due to infrared 
radiation rather than collision for our T v ;b and T rot . This 
is because collisional excitation would make T rot at V2—I 
equal to T v ;b. If the excitation were due to H2 collision 
at ^200 K the gas needs to be at or above the critical 
density for HCN y 2=1 , which is nn 2 ~ 5 x 10 11 cm~ 3 
(jZiurvs fc Turnerlll986t ). Radiative excitation is energet- 
ically possible. Each rotational state J at the vibrational 
ground level has associated infrared absorption lines in 
the P, Q, and R branches. The absorption lines supply 
energy to the HCN(J,J— 1) emission line, and weaker ro- 
tational lines at «2=1, in the limit of radiation-dominated 
excitation. Thus, an energy requirement for IR-pumping 
is 

^ fi Sir/Air > S mt /X TO t, (1) 

i=P,Q,R 

assuming that all the lines have the same velocity width. 
Here, Sir is the continuum flux density around the ab- 
sorption wavelength Air, fi is the fractional absorption 
depth, and S ro t is the mean flux density of the rota- 
tional emission line at the wavelength A ro t- Parameters 
for our case are Atr =14 jj,m, S'i4 Aim =2.34 Jy and Jq ~ 0.3 
(jLahuis et al.l 12007). and A rot =850 ^m and S rot < 1 Jy 
for J=4. The condition of Eq. (JXJ) is satisfied with a wide 
margin; the excess energy heats gas through collision. In 
this radiative excitation, the observed T v ib indicates that 
the brightness temperature of the excitation source seen 
from the gas must be > 200 K at 14 /xm. 

The vibrational excitation significantly affects rota- 
tional excitation in the vibrational ground level. We see 
this in a model with two vibrational levels, v=0 and 
1, whose energy gap is To in temperature. Each vi- 
brational level has rotational sub-levels labeled with J. 
For a molecule at (i;,J)=(0,J) there exist approximately 
e -T /T vib mo i ecu i es a t v —\ that will spontaneously decay 
to the (0,J) state. Thus the rate of spontaneous transi- 
tion to (0,J) from v—1 is e~ T °/ Tvib A V jb, and radiative 
excitation of the J level at v=0 through v=l needs 

e -T /T v , Mv . b > ^ 3 ^ Tv . b > Tq / ln Ax^L t ( 2 ) 

/ ^rot J 

where A v ;b and A Iot are the Einstein A coefficients for the 
vibrational and rotational transitio ns, respectively. This 
is equ ivalent to Equation (6) of iCarroll fe"~ Goldsmith 
1(1981 but is more useful when T v jb is available. Sim- 
ilarly, the radiative excitation through v=l exceeds col- 
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lisional (de)excitation among v=0 rotational levels when 

e- T °/ T - ib A vih >n OJ i7i 1 j- ll (3) 

where n co \ is the number density of the main collision 
partner and 7j,j-i is the collisional rate coefficient. In 
other words, radiative pumping to v—1 is equivalent to 
a gas density e _T °/ Tvib ^4vib/7J,j-i in terms of the rota- 
tional excitation around J at v—0. Parameters for HCN 
are Th=1024 K and A rih=3.7 s _1 for the bending mode 
(jZiurvs fc Turnenll98"6h and 7j,j_i « 1 x 1CT 11 cm 3 s^ 1 
for HCN-H2 collisions at J<10 and gas k inetic temper- 
atures 30-300 K (jDumouchel et all 12010ft . For the ob- 
served T vib of 227 K, the rate e - T °/ T « il >A vih is 0.04 s" 1 
and 20 times larger (faster) than A TOt j = 4_3. The mini- 
mum T vib to satisfy Eq. © is 86, 137, and 218 K for 3=1, 
4, and 10, respectively. Thus the observed Tab satisfies 
Eq. ([2]) for Js up to about 10. The IR-pumping through 
i)2=l matches nn 2 ~ 10 9 cm" for our T v ib according 
to Eq. although this is somewhat overestimated be- 
cause collisional transitions are not limited to |AJ| = 1 
and half of the vibrational transitions are in the Q branch 
where A J = 0. 

The relatively low rotational temperature 1^(^2=1 ! 
J^4)=31±13 K probably suggests that the system is not 
in local thermodynamic equilibrium (LTE). A caution 
here is that T rot may be higher and even comparable 
to T v ib depending on the uncertain 355 GHz continuum 
mentioned in |j3j T rot would be 59ll\j 5 K if we adopt the 
lower continuum level estimated from the LSB. A possi- 
ble reason for the lower rotational temperature than T v ib 
is that rotational excitation at low J (< 5) is strongly af- 
fected by collision and (sub)millimeter radiation; note 
that HCN(4-3) is optically thick. The T rot at v 2 =l 
may be relatively low accordingly. Rotational temper- 
ature can be higher in higher J since collisional excita- 
tion in <100 K gas gives way there to IR-pumping as 
the main excitation mechanism. This may explain the 
HCN (rotational) temperature of 300 K (±30%) that 
lLahuis et al.l (|2007ft estimated from the Q-branch absorp- 
tion profile under LTE assumption, because high J pop- 
ulation broadens the profile. 

For comparison, we analyzed the v = HCN lines 
alone using the large- velocity-gradient mode l via RADEX 
(jvan der Tak et al.ll2007t iSchoier et al.l 12005ft and found 
two types of solutions. We used for this the abun- 
dance ratios m e ntione d above and HCN (1-0) data from 
llmanishi et alJ (|2004l ). One solution has thermalized 
HCN at 3=4, N HC n/AV 10 14 7 c m - 2 (km s" 1 )" 1 , 
Tkin~40 K, and nn 2 > 10 7 5 cm~ 3 . Another has subthcr- 
mal excitation at J=4, moderate opacity of H 13 CN, and 
such parameters as jVhcn/AV « 10 15 cm _2 (km s -1 ) -1 , 
Tkin> 30 K, and nn 2 ~ 10 6 cm~ 3 . These temperatures 
are consistent with our T rot at v%=l and the high column 
densities are in accord with the deep obscuration toward 
the nucleus. However, these solutions cannot explain the 
observed HCN V2—I lines. 

To summarize, the vibrationally excited HCN is most 
likely IR pumped, and the observed degree of vibrational 
excitation suggests strong influence of the excitation path 
through «2 = 1 to the rotational excitation in the ground 
state although it may not dominate low-J excitation. 
Properties of (dense) molecular gas in NGC 4418 cannot 



be studied through HCN ignoring the IR-pumping. In 
further study, one needs to reexamine our simplistic as- 
sumptions that the emitting regions of «2=0 and 1 lines 
are cospatial and have a single characteristic tempera- 
ture and density. In reality, gas density must be diverse, 
and there must be low density gas where HCN is excited 
only under the IR-pumping. Thus the total flux of any 
HCN v=0 line should be enhanced at least for that com- 
ponent. Another improvement in modeling would be to 
includ e atomic hydrogen and electrons as c ollision part- 
ners (|Scoville et al.lll980t iFaure et al.ll2007ft . 

5. OTHER LINES AND OTHER GALAXIES 

HC3N is detected in its vibrational ground and excited 
states though some transitions are marginal or blended 
with other lines. Its population diagram in Fig. [4] shows 
that the measured line fluxe s are consistent with those o f 
other transitions detected by lCostagliola fc Aaltol (|2010ft . 
The diagram indicates overall excitation temperatures in 
the range of ~1 00-300 K, except fo r the E /k < 100 
K levels where iCostagliola fc Aaltol (|2010ft estimated 
T rot =29 K for v=0. They also obtained T vib ~500 K. 
These temperatures are compatible with our HCN exci- 
tation temperatures. 

The HCN to HCO + ratio of the mean brightness tem- 
peratures is 1.65 ± 0.07 for J=4-3 and 1.61 ± 0.08 for 
J=3-2 on the assumpti on of the same em itting vol- 
ume; it is 1.9 for J=l-0 (jlmanishi et al.ll2004D . Brighter 
HCN than HCO + in some Seyfert galaxies and in this 
galaxy has been proposed to be a result of or evidence 
for ch emistry driven by X-ray radiation from active nu- 
cleus (jKohno et al.l 1200 It llmanishi et ail 12004ft alth ough 
there have been d ifferent views (e.g.. lLahuis et~alll2007l: 
iBaan et al.l 12010ft . It is conceivable at least for objects 
as deeply embedded as the nucleus of NGC 4418 that 
the IR-pumping also affects this ratio. The minimum 
brightness temperature of the infrared source needed for 
the radiative pumping of J=l is given by Eq. ^ as 
T b mm = To/ ln^vib/A-ot 1-0) because vibrational tem- 
perature is at or below the source brightness tempera- 
ture at the wavelength of the vibrational transition. This 
minimum temperature is 86 K at 14 /im for HCN and 
156, 115, 107, 80, and 33 K at 4.6, 7.8, 12, 15, and 
45 ^m for CO, CS, HCO+ , N 2 H+, and HC 3 N, respec- 
tively (IChandra et al.|[T996l 119951: iMauclaire et~aL1ll995l: 
iHeninger et all 120031 . for A vih ). We note that HCO+ 
needs higher Tb than HCN for IR-pumping and that vi- 
brationally excited HCO + (3-2) is not detected at an ex- 
pected frequency of 266.8 GHz (see Fig. IR-pumping 
may therefore play a role in NGC 4418 for the high 
HCN-to-HCO+ ratio and maybe also for HCN-to-CO ra- 
tio and the bright N 2 H+. Its synergy wit h chemistry in 
hot galaxy nuclei fe.g.. lHarada et aU2010ft needs further 
study. 

We searched published spectra for the HCN(t>2 = l 1 ^) 
rotational lines at +400 km s^ 1 of HCO + (w = 0) to as- 
sess the prevalence of the HCN infrared pumping. We 
found no convincing case though a possi ble sign is at 
NGC 6240 nucleus in IWilson et ail (|200g ). The limita- 
tions are sensitivity, spectral baseline, and blending with 
the brighter HCO + lines. The lines must be as narrow 
as in NGC 4418 for unambiguous detection. 

6. CONCLUDING REMARKS 
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We evaluated the effect of infrared pumping to the ro- 
tational excitation of HCN using vibrational and rota- 
tional temperatures. The vibrational temperature be- 
came measurable with our detection of rotational lines 
from vibrationally excited HCN. The high T vib in NGC 
4418 suggests that HCN is IR pumped and cannot simply 
trace dense gas in the nucleus. Other molecules should 
also be subject to the pumping in different degrees ac- 
cording to their energy structures, A coefficients, spatial 
distributions, and the spectrum of the radiation field. 
Observations of T v ib in various molecules and galaxies 
will tell us the overall significance of IR-pumping in ex- 
tragalactic observations. Regarding NGC 4418, our high 
T v ib of HCN agrees with earlier studies to suggest a > 200 
K excitation source in the nucleus. The (sub) millimeter 
lines from vibrationally excited HCN are our new tool to 
probe hot molecular gas close to the nuclear heat source. 

We thank Daniel Espada for the CO(2-l) spectrum, 
Mark Gurwell for advice on flux calibration, the SMA 
operation team for the service observations, and the ref- 
eree for helpful comments. This research extensively used 
the NASA/IPAC Extragalactic Database, NASA's As- 
trophysics Data System, the Cologne Database of Molec- 
ular Spectroscopy, and the Splatalogue database. This 
work was supported by the grant 99-2112-M-001-011- 
MY3 from the National Science Council of Taiwan. 
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Table 1 

Emission Properties of NGC 4418 



Emission 
(1) 


f t 
[GHz] 
(21 

1-7 


E u /k 
PK1 


Ffline") S\ 

O V JVIIl a 

(4) 


°th ""tot 

— , mjy] 
(5) (6) 


note 
(7) 


HCO+(J=4-3) 


356.734 


42.8 


82 


3 


9 




HCN(t) 2 =l lf ,J=4-3) 


356.256 


1067.1 


30 


3 


1 




HC 3 N(>7=l lc ,J=39-38) 


355.566 


662.2 


6 


3 


3 


(a) 


HC 3 N(J=39-38) 


354.697 


340.5 


25 




10 


(b) 


HCN(J=4-3) 


354.505 


42.5 


133 


3 


11 


(c) 


CO(J=3-2) 


345.786 


33.2 


746 


3 


75 


(<i) 


CO(J=3-2), 10" 


345.796 


33.2 


782 


3 


78 


(d) 


H 13 CN(J=4-3) 


345.340 


41.4 


41 


1 
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Note. — (2) Line frequency. (3) Upper level energy, 
km s" 1 and continuum flux density in mjy in the central 6" 


(4) Line flux in Jy 
(1 kpc). Line fluxes 



are integrated over 1900-2300 km s _1 unless otherwise noted and are continuum 
subtracted. No correction is made for missing flux. CO (3-2) flux is also given for 
the central 10". No other emission was resolved by the 6" beam. (5) Thermal 
error of (4). For lines, this includes thermal noise of continuum propagated 
through continuum subtraction. (6) Total error of (4) including <r t h and 10% 
error of flux calibration, (a) Marginal detection with possible blending with 
weaker lines, (b) Large uncertainty due to line blending, (c) No subtraction 
is made for the blended HC 3 N(39-38). (d) There is likely contam ination by 
HC 3 N(J=38-37) at / rcst =345.609 GHz (Costagliola & Aa'lto 2010). (e) This 
line is blended with CO(3-2) and likely with HC 3 N(38-37). Its flux is obtained 
by doubling the line flux in 2100-2300 km s _1 assuming a symmetric profile, 
(f) Flux from 1900-2197 km s . (g) Continuum flux densities are measured 
outside the 1900-2300 km s -1 velocity ranges of the lines marked in Fig. [2] 
Also excluded for 0.85 mm continuum are HCN(4-3) velocities of 2300-2320 
km s — 1 . (h) Flux density is measured at the HCO + (3-2) velocities below 1900 
km s — 1 . 
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Figure 1. HCN energy level diagram, with i-doublets magnified 10 times on the right. Blue arrows are the rotational transitions that we 
detected. Allowed vibrational transitions are i>2 = <-> l 1 * , AJ = (Q branch) and V2 = <-¥ l le , AJ = ±1 (P and R branches); some 
of them are plotted. 
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Figure 2. Spectra of the central 6" (1 kpc) of NGC 4418. Major lines are marked at the redshift of the galaxy (2100 km s ). Full 
line names are in Table [l] Velocity axes are given for strong lines. Baselines of the shaded parts of the spectra are the continuum levels. 
Pairs of arrows in the top-right and bottom-left panels mark the frequencies of the HCO + (t)2=l) doublet. Scaled spectra of 3C273 from 
the same observations are also plotted to verify passband calibration. 
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Figure 3. Comparison of CO(3-2) and CO(2-l) spectra, showing that CO(2-l) has no emission (i.e., no high-velocity gas) at 2500 
(=V S y S + 400 ) km s —1 . The con tinuum-subtracted spectra of CO(3-2) and (2-1) are from the data in Fig. [2] and 3" resolution SMA 
observations (Espada ct al. 2010), respectively. 
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Figure 4. Population diagram of HC3N. The lower and upper abscissae are rotational energy measured from the rotational ground level 
of each vibrational state and rotational quantum numbe r J, respectively. Th e ordin ate is population divided by the degree of degeneracy, 
Nj v i/gj, in arbitrary unit. Open data points are from Costagliola & Aalto (2010) and filled ones are from our observations. Error bars 
are ±la. No correction is made for optical depth, but shifts needed for the correction are shown in the right side of the plot. The data 
for vj = 1 lines (red) are shifted against those for the vibrational ground ones (blue) using a vibrational temperature T v i^ = 300 K. (The 
V7=l state is 321 K above v=0.) Shifts for other T v ;ij are shown in the top-right corner. To its left are the slopes of rotational population 
for various rotational temperatures. 



